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Abstract 



Tensor models are more-index generalizations of the so-called matrix models, 
and provide models of quantum gravity with the idea that spaces and general 



relativity are emergent phenomena. In this paper, a renormalization procedure 
for the tensor models whose dynamical variable is a totally symmetric real three- 
tensor is discussed. It is proven that configurations with certain Gaussian forms 
are the attractors of the three-tensor under the renormalization procedure. Since 
these Gaussian configurations are parameterized by a scalar and a symmetric two- 
tensor, it is argued that, in general situations, the infrared dynamics of the tensor 
models should be described by scalar-tensor theories of gravity. 
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1 Introduction 

In the scenario of emergent spacetime [H El El IU EJ EJ [7] , the notions of spacetime and distance 
are not fundamental, but they are just infrared phenomena which appear only after underlying 
dynamics of fundamental degrees of freedom are averaged over. Therefore, to compare the 
infrared dynamics of such a scenario with the conventional spacetime formulation of physics 
such as the general relativity, it is necessary to give an appropriate definition of distance 
and derive the metric tensor field on an emergent spacetime. An explicit example of such a 
definition is given in the noncommutative geometry, where the notion of distance is defined 
through a Laplacian or a Dirac operator on a noncommutative spacetime [8]. Another example 
is given in the simplicial quantum gravity, where such operational definitions of distance have 
been used for comparisons with the general relativity [7J. 

Because of the infrared character of the notion of distance, different operational definitions 
of distance are not generally identical in short distances but converge in long distances. This 
implies that a metric tensor field on an emergent spacetime could be derived more directly 
from application of a renormalization procedure to underlying small distance dynamics of 
fundamental degrees of freedomo In this paper, I will discuss a renormalization procedure 
for the tensor models, which are models of emergent space, and see how a metric tensor field 
emerges. 

The tensor models were originally proposed to describe the simplicial quantum gravity in 
more than two dimensions [TUl [TT] . Although the subsequent developments [121 021 El HSl 
[TBI nT] have demonstrated some interesting connections to the loop and spin foam quantum 
gravities, the tensor models seem to have not yet become a fully successful approach under the 
original interpretation of the tensor models based on the diagrammatic relations between the 
tensor models and simplicial manifolds in more than two dimensions. In view of this status, 
the author has proposed another interpretation of the tensor models as theory describing 
dynamics of fuzzy spaces [18]. This reinterpretation has the following advantages compared 
with the original one. Firstly, a three-tensor is enough to describe all the dimensional spaces. 
Secondly, the semiclassical analysis of the tensor models can be considered to have the physical 
meaning that a classical solution and small fluctuations around it correspond to a background 
space and field theories, respectively. Lastly, the diffeomorphism invariance of the general 
relativity can naturally be incorporated into the symmetry of the tensor models. Indeed, the 
tensor models which are fine-tuned to possess classical solutions with certain Gaussian forms 



*In fact, in the context of the AdS/CFT correspondence in string theory, coarse-graining and renormal- 
ization are argued to be essentially important in relating boundary field theories to bulk classical geometries 
0. 



have numerically been analyzed, and the agreement between the tensor models and the general 
relativitjiiJ has been obtained in various dimensions [T9l 



The purpose of this paper is to give a first step to generalize the above results of the 
fine-tuned tensor models to more general tensor models. It will be proven that the Gaussian 
configurations are the attractors of the three-tensor under a renormalization procedure. Since 
these Gaussian configurations can be parameterized by a scalar and a symmetric two-tensor, 
it will be argued that, in general situations, the infrared dynamics of the tensor models will 
be described in terms of scalar-tensor theories of gravity. 

Coupled field theories of the metric tensor and a scalar have been considered for various 
physical reasons [24]. For example, motivated by Mach's idea on inertial induction, the Brans- 
Dicke theory [25] has a scalar field which plays as a variable gravitational "constant" . In string 
theory, the dilaton field is necessarily introduced to keep the worldsheet conformal invariance 
after quantization [26]. Scalar fields are actively discussed to model the dark energy and the 
dark matters, which have clearly been shown to exist in the universe by recent astrophysical 
observations [27] . Therefore it will be highly interesting, if it is concluded that the tensor 
models necessarily generate a scalar field as well as the metric tensor field in their infrared 
dynamics. 

This paper is organized as follows. In Section [2] the tensor models in this paper and the 
Gaussian configurations are introduced. In Section [3], a renormalization procedure for the ten- 
sor models is introduced. In Section [4] it is proven that the Gaussian configurations are the 
attractors of the three-tensor under the renormalization procedure. In Section [5] two numer- 
ical examples are given for showing how configurations converge to Gaussian configurations 
by iterations of the renormalization procedure. In Section [HI it is argued that the infrared 
dynamics of the tensor models should be described in terms of scalar-tensor theories of gravity. 
The last section is devoted to summary, discussions and future directions. In Appendix [A] 
some details of the proof in Section [4J are provided. 

2 Tensor models and Gaussian configurations 

The general idea of regarding the tensor models as theory describing dynamics of fuzzy spaces 
is based on the observation that a fuzzy space is characterized by an algebra of functions on 



' Here the terminology "general relativity" is not appropriate in its strict sense. I use it to mean a theory 
which has a metric tensor field with Euclidean signature as its dynamical variable, and has the diffeomorphism 
symmetry. Moreover, its action may not be the Einstein- Hilbert action, R, but may be R , or other invariants. 



it, f a fb = C a b c f c [IBjIil By considering the three-tensor C to be a dynamical variable, one can 
obtain the tensor models as theory of dynamics of fuzzy spaces. Depending on the detailed 
choices of a dynamical variable and a symmetry, there exist various kinds of tensor models. 

In this paper, the dynamical variable of the tensor models is a totally symmetric real 
three-tensor C a bc (a,b,c — 1,2, ... ,N), which satisfies 

Cabc ^ K, 

1 

^abc ^bca ^cab ^bac ^acb ^cba- 

The tensor models are assumed to be invariant under the orthogonal transformation, 

C' abc = M a a 'M b b 'M/C a , b , c ,, M e O(N). (2) 

Then an example of an invariant action can be given by 

o y® /~i /~i i yi /~i /~i /~i /~i y^ /~i /~i /~i ^ fo\ 

iJ — — ^abc^abc "T — ^abc^dbc^def^aef T~ ^abc^ade^bdf^cef , \<JJ 



where the repeated indices are assumed to be summed ovei[f|, and g^ are real coupling constants. 
This is one of the simplest actions, which is known to have a number of classical solutions 
corresponding to various dimensional fuzzy spheres and tori [28]. In Section [51 I will use 
its classical solutions as examples for applying the renormalization procedure introduced in 
Section [3J 

Next, let me introduce the Gaussian configurations [TH [201 [21] . m concrete applications, 
the indices of the three-tensor are taken to be discrete and finite, but, as ideal cases, let me 
assume the indices to be continuous and infinite in the following discussions. In the coordinate 
basis, the Gaussian configurations are defined to take the form, 

C 9 X1X2X3 = 5exp [-0 ((a* - x 2 f + (x 2 - x 3 ) 2 + (s 8 " *i) 2 )] , (4) 



where Xi are .D- dimensional coordinates Xi = (x},x 2 , . . . ,xf), B and (3 are positive real 
constants, and x 2 = g liV x ll x v . Here g^ v is a constant real symmetric two-tensor and is assumed 
to be positive-definite to demand exponential damping of C 9 at large mutual distances among 
Xi. The parameter /3 is redundant in the sense that it can be absorbed into the rescaling of 
g^v, but it will be kept for the later convenience to characterize the order of the fuzziness, 
1/vft while gp,~0(l). 



-^The algebra considered is generally nonassociative. Otherwise, an algebra cannot encode geometry of a 
fuzzy space. For example, geometry is encoded in a Laplacian (or a Dirac operator) on a noncommutativc 
fuzzy space 8 , but not in a noncommutative algebra itself. 

§This will be assumed in the rest of this paper, unless otherwise stated. 



Because of the obvious translational invariance of (jlj), it is more convenient to express (jl]) 
in the momentum basiqj. By applying Fourier transformation to the indices Xi of (jl]), one 
obtains 

C* ip2P3 = Aexp [-a {( Pl ) 2 + (p 2 ) 2 + (p 3 ) 2 )} S D ( Pl + p 2 + p 3 ) 7 (5) 

where a = 1/(12/3), A is a positive constant proportional to B, pi = (pj,p 2 , . . . ,pf), and 
p 2 = g^PuPu with g^ v being the inverse of the matrix g^. 

The algebra of functions on a usual continuous space is given by f xl f X2 = 8 D {x\ — x 2 ) f xi . 
This corresponds to C XlX2X3 = 5 D (xi — x 2 )5 D (x 2 — x 3 ) by identifying C a bc with the structure 
constant of the algebra, f a fb = C a t, c f c . In fll]), the delta functions are smoothened by Gaussian 
functions. In addition, (J3J) (or (JSJ) has the Poincare invariance. Thus the three-tensor (j4]) (or 
(E}) can be considered to represent a D-dimensional fuzzy infinite flat space [29] . 

3 Renormalization procedure 

The philosophy of the Kadanoff- Wilson renormalization procedures is in iterations of coarse 
graining processes. In a discrete system, such a process is generally a discrete step of defining 
new renormalized dynamical variables by averaging over contributions of underlying dynamical 
variables on nearby sites. On a fixed lattice, the meaning of nearby sites can obviously be 
defined by giving a renormalization scale, but in a gravitational system, one has to take into 
account the fact that a scale is a dynamical quantity determined by dynamical variables. 
Hence, in analogy with a gravitational system, a renormalization step in the tensor models 
should be expressed solely by the dynamical variable, and it will generally take a form, 

Cl = R(C) abc , (6) 

where R(C) a b c is a totally symmetric three-tensor as a function of C, and C R is a new renor- 
malized dynamical variable obtained after one step of a renormalization procedure. The O(N) 
symmetry (J2J) is assumed to be respected in (jSJ). 

In fact, there exist various possibilities of a renormalization procedure R(C), but I assume 
that the qualitative behavior in the infrared limit will not depend on the details of the choice 
of R(C), if it is taken appropriately. The simplest choice of R(C) [20] is as follows. Let me 
consider a configuration of C which is not far from the Gaussian form (J3J). Then, if x,y are 
nearby points, C xyz f z gives a function averaged over within the distance of order 1/a/5 from 
x, y. Similarly, one can consider another averaged function C x > y ' z 'f z ' 7 and take the product 



^In the momentum basis, the contractions of indices must be taken for pairs of p and —p. 



of these averaged functions, C xyz f z C x / y t z 'f z i = C xyz C x ' y > z 'C zz ' w f w , which will determine a 
renormalized structure constant. Contracting the unwanted indices y and y', one obtains a 
renormalized three-tensor as 

R{C) a bc = C a d e CbdfC ce f. (7) 

Note that this expression respects the O(N) symmetry. 

It is obvious that there exist various possibilities to insert C a b c into ([7]) to change the 
detailed way of averaging over. I assume that these changes do not alter the essential properties 
of the infrared behavior of a renormalization procedure. I will show a numerical evidence for 
this robustness in Section In the following section, I will discuss the simplest renormalization 
procedure ([7]). 

In general, the renormalization procedure ([6]) does not change the number of degrees of 
freedom, and hence this is not the same as the usual Kadanoff- Wilson renormalization proce- 
dure. But as will be discussed in Section [6], after a number of iterations of the renormalization 
procedure, the degrees of freedom can be accumulated into a scalar and a two-tensor field. 

4 Proof of convergence to Gaussian form 

Let me first discuss the uniqueness of the Gaussian configurations ([5]) up to the second order 

of momenta in the exponential, under the assumption of the momentum conservation. Since 

C PlP2P3 must be totally symmetric, the most general expression up to the second order is given 

by 

C PiP2P 3 =Aexp[h IJ, pi ll + h^p 2fl + h^p 3fl 

+ h^p l ^p 2v + h^p 2 ^p 3u + h^p 3 ^p lv 

+ h^ Pl , Pl „ + W v p 2 ,p 2 „ + W v p 3 M p 3 u] 

x 5 D (pi+p 2 +p 3 ), 

where W ,v and W are symmetric, and the momentum conservation is expressed by the delta 
function. Because of p\ + p 2 + p 3 = 0, the sum in the first line vanishes. Moreover, by using 
h^ v {p\ +p 2 +P3)fj,(pi +p 2 +Pz)u = 0, the expression in the second line can be converted to that 
in the third line. Thus ([5]) is the most general expression up to the second order of momenta 
in the exponential. 

Now let me apply the renormalization procedure (|7j) to the Gaussian configuration (JSJ). 
Noting that the contraction of indices must be taken for pairs of p and — p in the momentum 



basis, one obtains 

•tt(P ) PlP 2P3 = 

= A 3 5 D ( Pl +p 2 + p 3 ) / d D q exp (-a(p? + p 2 + #| + 2 g 2 + 2(p x + g) 2 + 2(p 2 - g) 2 )) 



/^g (-19 n9 

Pi 9i 92 VI -91 93 PZ -Q2 -93 



A s 



* D g 

D 



(6a; 



exp ( --a(pl +p\+ p\) ) S u (p! +p 2 + pa) 



(9) 
where g = Det(^i/)U!I This expression is again a Gaussian configuration, and the parameters 
A a are renormalized as 



)fi 



A fl = A» 






(6a 



(10) 



«b = -a. 
3 

This proves that the Gaussian configurations ( 

malization procedure (jTJ). 



or (jSJ) are on the trajectories of the renor- 



In the following, I will show that the Gaussian configurations are the attractors of the 
three-tensor under the renormalization procedure. For simplicity, let me assume that g^ iv = 
9^ v = $nu- This simplification does not ruin the generality of the discussions below, because a 
linear transformation of the momenta in (jSJ) can transform a general positive-definite g^ v to 
5^. Let me consider a small deviation from the Gaussian configurations (jSJ) in the form, 

C vip2pz = A ( l + f (V®PuVaP2,Vap 3 ))exp(-a(p 2 1 +pl+p 2 3 ))5 D (p 1 +p2 + P3), (11) 

where the function / represents a deviation from (J5J), and must be symmetric with respect to 
its arguments. Note that this deviation does not violate the translational invariance or the 
momentum conservation. The violating case will be discussed in Section EJ 

Now let me discuss the renormalization of the deviation /. Let me assume that the de- 
viation is so small that the renormalization procedure can be approximated by the linear 
approximation in /. The linear approximation is enough to see whether the Gaussian con- 
figurations are attractive or not. Applying the renormalization procedure (jTJ) to (ITT)) , one 



II Note that g is for g^ u , while p 2 — g^p^Pv 



(12) 



obtains 

-ti{^ ') 'pip2P3 — ^Pl<7l q'2 < ~ y P2 — 91 93 ^"P3 -<?2 -<J3 

~ A 3 5 D ( Pl + p 2 + p 3 ) / cPg [1 + / (Vapi, V^q, y/a(-q - Pi)) 

+/ (v 7 "^, -v 7 "?, v 7 "^ - P2)) + / (Va?3, Va(p 2 - q), Va(q + pi))] 
x exp (-a(p? + P 2 + P3 + V + 2 (pi + <?) 2 + 2(p 2 - g) 2 )) 

= A 3 exp I --a(p\ +pl +pl) ) 5 D (pi + p 2 + P3) / d D gexp(-6ag 2 ) 

x [l + / (y/api, y/a(q + (p 2 - Pi)/3), \/a(-? + (P3 - Pi)/3)) 
+ / (a/«P2, \/a(? + (P3 - P2)/3), \/a(-9 + (Pi - P2)/3)) 
+/ (\/aP3, v 7 "^ + (Pi - P3)/3), ^fa(-q + (p 2 - P3)/3))] , 

where only terms linear in / are taken from the first to the second equalities. After replacing 
A, a with the renormalized Ar, an in ( ITUj) . the renormalized deviation Jr can be read from 
(H2D as 

t ( 1— I— /— \ [ 18a A^ f ab ( 18a fl 2 \ 
Jr [VoiRPi, VorP2, VOrP3) = I -z — J / d q exp I — q I 



+ (cyclic permutations of pi,p 2 ,p 3 ) 



In this paper, I consider an / which contains only a finite number of polynomials of pi. In 
the coordinate basis, such a deviation corresponds to adding a finite number of polynomials 
of Xi — Xj to the constant B in (J3J). Therefore such a deviation is local in the sense that it 
does not change the exponential damping behavior of C XlX2X:i at large mutual distances among 
Xi. A general deviation with such damping behavior may have an expression with an infinite 
number of polynomials. But, because of the damping behavior, it will be enough to consider 
only the central region, and it will be possible to approximate a deviation with good accuracy 
by a finite number of polynomials. Then the proof below will be applicable also for such 
general cases. 

Because of the momentum conservation pi + p 2 + p 3 = 0, the expression of / in terms 
of pi is not unique, that causes difficulties in the proof. This ambiguity can be avoided by 
expressing / solely in terms of two of the momenta. Since / is symmetric with respect to the 
momentum arguments, there exist three equivalent expressions, each of which can be obtained 



from each choice of two momenta, as 



f(y/api,y/ap 2 ,y/ap 3 



Pl+P2+P3=0 



m+n<M 



f(y/api, v^Pa, y/a(~Pi ~ Pa)) = Yl «^/i" 1 ^"" ,M ^>i (ll • 

m,n=0 

m + n^l 

f(V^(-p 2 - p 3 ), V«P2, v^Ps) = ^ o: s ^V 1 -'^ ;i/1 "' / "p 2/ , 1 ■ 



m,n=0 
m + n^l 



•Plu m P2vi • " "Pa* 



■P2 i i m PZv 1 



■PlVn 



(14) 



where M is the maximal degree of the polynomials of pj, and /i are the coefficients. These /i are 
symmetric with respect to the indices: /^i-MmM-^ — } l yi-v n ;nx...nm an( j ^(lj-^wiVtii-Vw — 
j l ^i...^i m ;ui...v n ^ ^gj-g ^ a ' denote permutations. Here the right-hand sides of f|T4|) do not contain 
any linear terms with m + n — 1, because of the momentum conservation, Pi + P2 + P3 = 0. 

Substituting the polynomial form f[T4"l) into (fT3l) . the integration over g of each term will 
generate terms of lower or equal degrees. Therefore the renormalization procedure closes with 
the form in (TT4"|) . and it can be expressed as a linear transformation with an upper triangular 
matrix, 

/ h-Ro \ I Roo R02 
R22 






h 



R2 



\hRMj 



Rom \ I h \ 
h 2 



R2M 



\ 



Rmm) 



(15) 



\huj 



where hi and hm symbolically represent the coefficients of the terms of degree i in / as in (1141) 
and fa, respectively, and Rij symbolically represents the sub-matrix relating hj to hm- The 
eigenvalues of the linear transformation (TT5|) determine whether the Gaussian configurations 
are attractive or not under the renormalization procedure. Because of the triangular form 
of the linear transformation, the eigenvalues are determined solely by the eigenvalues of the 
sub-matrices Ra on the diagonal. 

The matrices Ra on the diagonal represent the contributions from / to /r which does 
not change the degrees of the polynomials. These contributions can be evaluated simply by 
ignoring q in the arguments of / in (113I) . and is explicitly given by 




3orP 3 -Pi \ , v ... , , ( 16 ) 

+ (cyclic permutations of Pi,P2,P3J- 



To discuss the matrix Rkk [K < M), let me assume that / (and f e ^ ua ) contain only polyno- 
mials of degree K. Then ffT6|) becomes 

K_ 

jequai (j) liV2iV ^ = f - J f (pi, 2 - , 3 J + (cyclic permutations of Pi,p 2 ,P 3 ). 

(17) 

To discuss the qualitative behavior of the process (TT7|) . let me introduce a norm of / defined 
by 

H-f K «, M £7 «J/(Pi>ft»>**)l> ( 18 ) 

(P1) 2 >(P2) 2 ,(P3)^<1 
Pl+P2+P3=0 

where | | denotes an absolute value, and Max represents taking the maximum value under the 
conditions below it. From (TT7I) . by using |/i + /2I < |/i| + I/2I, one obtains 



j equal II < 3 ( £) " Max 

K 5/ (P1) 2 ,( P2 ) 2 ,(P3) 2 <1 

Pl+P2+P3=0 

< 3 I'l)' 11/11. 



, , P2- Pi P3 ~ Pi 
J [Pi, o > ^ 



(19) 



since the three momenta in the arguments of / in the first line of ( fl9l) satisfy 

. P2-P1 . P3-P1 n 
P ^^^ = 

/ ^2 (P2-P1Y (P3-Pl\ 2 ^i 
(Pi) , [~^j , [~^j < 1, 

under (pi) 2 < 1 and p\ +P2+P3 = 0. The inequality (fl9l) proves that, for i^ > 5, the deviation 
asymptotically vanishes by iterations of the renormalization procedure. This implies that the 
eigenvalues of Rkk {K > 5) are smaller than 1. 

For i^ = 0, 2, 3, 4, the inequality ( TT91) is useless because the numerical factor is larger than 
one. In these cases, one can explicitly write down all the possible forms of /, and explicitly 
apply ffTTj) to obtain the eigenvalues of the matrices Rkk- For K = 0, the only form of / is 
a constant and the eigenvalue is obviously 3. For K = 2, there exists only one kind of form, 
h^ v {piuPiu +P2uP2u +P3vP3u), and the eigenvalues turn out to be 1. For K = 3, a detailed 
analysis shows that there exists only one kind of form, and the eigenvalues are (3/5) 3//2 . For 
K = 4, there exist three kinds of form, and the eigenvalues are 11/25 and 3/25. The details 
of these computations are given in Appendix [A] 

The above discussions show that the eigenvalues of Rkk are less than one for K > 3, 
but not for K = 0, 2. This implies that, by iterations of the renormalization procedure, 



the fluctuations of order K > 3 asymptotically vanishes, while those with K = 0,2 do not. 
However, these remaining terms have the forms of a constant and h tJ,u (pi tl pi u + p 2il p 2v + 
P3^P3u), respectively, which can be absorbed into the redefinition of A and g^ v of the Gaussian 
configuration (jSJ), so that they are deleted from fluctuations. This proves that the Gaussian 
configurations (jl|) or (151) are attractive under the renormalization procedure (171). 

5 Numerical examples 

In this section, I will show two numerical examples of how configurations converge to the 
Gaussian oneso To test the renormalization procedure in concrete cases, the starting config- 
urations are taken from the numerical solutions to the equation of motion derived from the 
action ([3]). The equation of motion is given by 

- 00 Cabc + y {CadeCdefCf hc + (cyclic permutations of a, b, c)) - g 2 C ade C bdf C cef = 0. (21) 

This equation has a number of numerical solutions corresponding to various dimensional fuzzy 
tori and spheres [28J. In this section, as the simplest example, I will consider a solution corre- 
sponding to a one-dimensional fuzzy torus (a ring), which has a one-dimensional translational 
symmetry. 

Because of the translational symmetry, it is convenient to describe the solutions in the 
momentum basidlll. To consider a ring, the momentum index is assumed to take only in- 
teger values. A cut-off L is introduced so that the range of the momentum index is finite, 
p = —L, —L + 1, . . . , L. The one-dimensional translational symmetry can be realized by the 
momentum conservation of the solution, i.e. C PlP2P3 ^ only if p\ + p 2 + p% = 0. In the 
left figure of Figure El a solution for g = l,g 1 = 2,g 2 = 1, L = 100 is plotted. The form 
of the solution is very different from a Gaussian function. Nonetheless, after 10 iterations of 
the renormalization procedure, the configuration approaches to a Gaussian configuration as 
shown in the right figure of Figure [2j 

It will be interesting to check the robustness of the convergence to the Gaussian configu- 
rations under the changes of the renormalization procedure. As an example, let me consider 



"The Mathematica program used in this section can be downloaded from 



http://www2.yukawa.kyoto-u.ac.jp/^sasakura/codes/renormalization.nb. 

ft At first sight it is not evident whether the index p below has the physical meaning of momentum, since a 
space is dynamically generated and there is no preferred frame because of the O(N) symmetry ^. In fact, 
physical meaning of an index can only be understood by checking the internal consistency of an interpretation 
of a solution as a whole. For example, one can see that, in the basis x dual to p, the numerical solution respects 
locality, i.e. C X1 X2X3 is negligible unless Xj are nearby locations as in Figure [TJ This fact ensures that the dual 
basis x can appropriately be interpreted as a coordinate with a small fuzziness. 

10 




Figure 1: The solution C 0x2X3 representing a fuzzy ring for L = 10 is plotted in the coordinate 
basis. The values are normalized by Coo o- The coordinate x dual to p is defined by the discrete 
Fourier transformation J2 _ L exp(2iiipx/(2L + 1)). 
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Figure 2: The left figure is the plot of the starting configuration which is a solution of (1211) 
for go — 1, gi — 2, g2 — 1, L — 100. The right figure is the plot of the configuration obtained 
after 10 iterations of the renormalization procedure (J7|). The vertical and horizontal axes are 
C PlP2P3 in the momentum basis and p = \fp\ + p\ +p|, respectively. The values of C P1P2P3 
are normalized by C*ooo- The solid line in the right figure is a fit with a Gaussian function, 
c exp(-cip 2 ). 
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Figure 3: The configuration obtained after 10 iterations of the other renormalization process 
( 122]) . The solid line is a Gaussian fit. 



another renormalization procedure defined by 

R{C)abc = C ade C bdf C ceg C fhi C hig + (cyclic permutations of a, b, c) 



(22) 



which can be obtained by inserting some C to (J7|). Since (1221) is a sum of non-symmetric 
terms, this renormalization procedure does not seem to have the Gaussian configurations as 
stable configurations. Nonetheless, after 10 times of iteration, the configuration converges 
very well to a Gaussian function as shown in Figure O This fact suggests that the arguments 
of this paper might be more widely applicable than to the specific renormalization procedure 
(J7J). Presently, I have no proof of convergence for such general cases, but this character of 
convergence might be related to the so-called central limiting theorem in statistics. 

6 Effective field theory 

The discussions so far in the previous sections have assumed the translational invariance of 
the fluctuations. In this section, the fluctuations are allowed to have a very small momentum 
(pi +P2 +P3) 2 "C 1/a, and are assumed to respect the locality, i.e. C XlX2X3 is negligible unless 
Xi are mutually located within the distance of order l/y/Jf. 



iitv, 

§ 



Let me take the coordinate basis, and assume that the whole space is covered by patches 
of finite sizes well larger than the scale of fuzzyness, l/V/3 ~ a/«, but well smaller than the 
wavelength, 1/a/(pi + Vi + P3) 2 ^ \A*. Then because of the assumption of the locality of C, 
the renormalization procedure ([7]) is practically a local process within each patch. Moreover, 



^The two conditions are mutually compatible. In the coordinate basis, the first condition implies that 
fluctuations of C xx + ax +b have wave lengths of order l/-\/(pi + P2 + P3) 2 ^ \foi as functions of x, while the 
second one implies that C xx+ax +b is negligible unless a 2 ,b 2 < 1/(3 ~ a. 
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C can be considered to have the translational invariance with good accuracy in each patch. 
Thus the proof in Section H] will be applicable to each patch, and a number of iterations 
of the renormalization procedure jTf on C will result in a Gaussian configuration on each 
patch. Then the collection of the Gaussian configurations of all the patches will define a 
Gaussian configuration depending on locations in the whole space. Such a generalized Gaussian 
configuration may be expressed as 

C X1X2X3 = e^^^ + ^^g(x 1 ^g(x 2 ^g(x 3 ^ex P {-(5 {d{x Xl x 2 f + d(x 2 ,x 3 ) 2 + d(x 3 ,x 1 ) 2 )) , 

(23) 
where <fi(x) is a scalar field, g(x) = Det^^rc)), and d(xi,Xj) is the geodesic distance between 
Xi and Xj on a background of a metric tensor field 57^,(2). The above expression f[2"31 is 
covariant under the general coordinate transformation in the sense that the contraction of 
indices, C XlX2X C XX3X4 = j d D x^Jg~{x)e 2 ^ x " > ■ ■ ■ , is invariant, provided that <p(x) and g^x) 
transform covariantly as a scalar and a two-tensor field, respectively. The main motivation 
for the specific form (123]) is that this invariance is necessary to identify a part of the O(N) 
symmetry of the tensor model with a fuzzy analogue of the diffeomorphism symmetry. Note 
that the covariant expression ( 1231) is not unique: g^x) may be rescaled by using 0(x), e.g. 
9nv(x) — > ^^9iiv{x)i and there may be corrections in higher orders of 1//3, e.g. g^ u {x) — > 
g^ v {x) + R flu (x)/f3, etc. 

In the renormalization procedure ([7]), the number of degrees of freedom do not change, 
that is out of the philosophy of the Kadanoff- Wilson renormalization procedure. However, 
the above discussions show that if the above assumptions are met, the configurations can 
be well approximated by a scalar and a two-tensor field after a number of iterations of the 
renormalization procedure. Then the infrared dynamics of the tensor models will be described 
by an effective field theory of a scalar and a two-tensor field, where the latter field will be 
identified with the metric tensor field of the general relativity because of the fuzzy analogue of 
the diffeomorphism symmetry. Indeed, the form ( 1231) is a generalization of the one presented 
in the previous papers [20], [2TJ, which contains only g^ v (x) and has been used in the numerical 
analysis of the tensor models which are fine-tuned to have the Gaussian configurations as 
classical solutions. The comparison between the light modes in the tensor models and the 
general relativity has shown some remarkable agreements [191 1201 [21] . 

On the other hand, the discussions in this section show the possibility that a scalar field 
also appears in the infrared dynamics of the tensor models. Thus, in general situations, the 
infrared dynamics of the tensor models will be described in terms of scalar-tensor theories 
of gravity. Then it becomes a question why a scalar field has not shown up in the previous 
analysis of the fine-tuned tensor models. This is probably because, while the fuzzy analogue of 
the diffeomorphism invariance coming from the symmetry of the tensor model constrains the 
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two-tensor field to appear as a gauge field, the dynamics of the scalar field strongly depends on 
the details of the tensor models, and its appearance as a light mode is not always guaranteed. 

7 Summary, discussions and future prospects 

In this paper, I have discussed a renormalization procedure for the tensor models which have 
a totally symmetric real three-tensor as their only dynamical variable. It has been proven 
that configurations of certain Gaussian forms are the attractors of the three-tensor under the 
renormalization procedure. Then since these Gaussian forms can be parameterized by a scalar 
and a two-tensor, it is argued that, in general situations, the infrared dynamics of the tensor 
models will be described by scalar-tensor theories of gravity. 

An important assumption in the argument is the locality of a background classical solution 
and the fluctuations around it. This is not an evident assumption in the tensor models, 
because the three-tensor contains non-local elements. The locality of a background solution 
can directly be checked by explicitly evaluating the three-tensor. On the other hand, the 
fluctuations around it are generally not local. Then what is really assumed is that the non- 
local fluctuations are not important in the infrared dynamics of the tensor models. This is a 
problem of dynamics, which could be studied in numerical manners. 

In this paper, the dynamical variable of the tensor models is a totally symmetric real three- 
tensor. A natural extension of this is to consider a three-tensor satisfying the generalized 
hermiticity condition, which was considered in the original proposals of the tensor models 
[HI QUI EH] • This turns out to be equivalent to consider two real three-tensors, one of which 
is totally symmetric and the other is anti-symmetric. Thus this extension will inevitably 
introduce an anti-symmetric two-tensor field by similar discussions as in this paper. Then the 
whole content of fields will agree with the massless fields in the bosonic string theory [2H] . 
It would be highly interesting to study the properties of this new field to see whether this 
agreement is just an accident or has more meaning. 

Another interesting direction is to apply the renormalization procedure to the actual anal- 
ysis of the tensor models. As explained at the end of the previous section, it is presently not 
clear why a scalar field has not shown up in the previous analysis of the tensor models which 
are fine-tuned to have Gaussian classical solutions [HI [201 [2T] . Probably this is because the 
scalar field can appear as a light field only on limited conditions. In this regard, the tensor 
model with the action (J3J) will be interesting to study, because it is known to have various 
physically interesting classical backgrounds without fine-tuning of the coupling constants [28] . 
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The model will also be interesting as an example to study the characteristics of the scalar- 
tensor theories describing the tensor models. The forms of scalar-tensor theories of gravity 
depend on the physical motivations |24|. Among them, the idea of Brans-Dicke [25] seems very 
interesting in the present context, because, in the scenario of emergent spaces and gravity, the 
gravitational constant will be a dynamical quantity which will be determined not only by a 
background spacetime but also by matter distributions. 

A The eigenvalues of Rrk f° r K = 2,3,4 

In this appendix, for K = 2,3,4, I will explicitly obtain all the independent forms of the 
functions f(pi,P2,P3)<, which are symmetric with respect to Pi and are defined only on the 
momentum conserved plane, p\ + p 2 + p% = 0. Then f e ^ ua is computed by (TT7|) to obtain the 
eigenvalues of Rrk- 

The tensor h in this appendix is assumed to have the symmetries appropriate for (J14j) . 
i.e. \iV-\---v-m\v\---vn _ f l vi...u n ;m...tim anc j ^/V(i) ■ ■ /V(m) ;*V(i) ■ ■ ■'V(n) = j^i—iMn\vx—v n where er er' denote 

permutations. 

A.l K = 2 

Apparently, for K = 2, there exist two functions which are symmetric with respect to Pi as 

/l(Pl,P2,Ps) = K (pi^Plu+P2^P2u+P3^P3u) , 
/2(Pl,P2,Ps) = h%' V (pi lx P2„+P2 IJr P3v+P3nPlv)- 

However these are not independent, because, by substituting p 3 = —p\ — p 2 into ff24l) . these 
functions turn out to have the same form as 

fl(pi,P2,~Pl -Pi) =2K U (p lfl p lu +p 2ll p2v+PlvP2v), 
f2(pi,P2,~Pl -P2) = -hg"(pitfH v +P2 ll P2u+PlitP2v)- 



Thus there exist only one independent form. Substituting f\ into ( 1171) . one finds 

K_ 
flR al (PuP2,-Pl -P2) = f gj g/l(Pl,P2,-pi -P2) = fl(Pl,P2,-Pl~P2)- (26) 

Thus the eigenvalues of R 22 are 1. 
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A. 2 K = 3 

For notational simplicity, let me first introduce some short-hand notations, 

p 3 = h^fp&vPp, 

2 ( 27 ) 

(Pi,P2) = h^' p p lp p lu p 2p , 

and so on. These are convenient, because one can compute without taking so much care of 
the tensors h. The following formulas hold as usual: 

(pi + P2) 3 = p 3 i + 3plp 2 + Zpipl + pi 

((Pi+P2) 2 ,Ps) = (p 2 i + 2piP2+phP3) = (Pi,Pa) + 2 (piP 2 ,P3) + (pliPs), 

where p\p 2 = h pup p lfl pi u p 2p , and so on. With these notations, there exist two apparently 
independent forms, 

h = p\+pI + pI = -Kp\p2 + pipI), (29) 

h = (Pi,p 2 ) + (permutations oipi,p 2 ,p 3 ) = (pt,P2) + 2(p 1 p 2 ,pi + p 2 ) + (p 2 ,Pi), (30) 

where I have used ps = —p± — p 2 . In fact, f\ and f 2 are not independent, because, by 
substituting h IJ,u ' p = h pvp into (!30|) . one can derive the same form as (|29|) . Thus it is enough 
to consider only / 2 , and, by substituting it into (TT71) . one obtains 

#'SN "2" 

requal < " 



Thus the eigenvalues of i? 33 are (3/5) 3 / 2 . 

A.3 K = 4 

With similar simplified notations as above, there exist three apparently different forms, 
/1 = p\ + P 2 + Pt = 2(Pi + 2p 3 p 2 + 3p 2 p 2 + 2 Pl ^ + pl) } 
/2 = (Pi,P2) + (permutations of pi,p 2 ,p 3 ) 

= -2(p?, Pl ) - 2(plp 2 ) - (plp 2 ) - (pl Pl ) - 3(p 1 p 2 (p 1 +p 2 ),p 1 +p 2 ), (32) 

/3 = (Pi^pI) + (cyclic permutations of Pi,P2,P3) 

= (pI,P 2 i) + (pIpD + 3(p 2 ,P 2 ) + 2(piP 2 ,Pi +pl), 

where I have used (Pi,p 2 ) = (p 2 jP?)> "which comes from h pu ' ,pcr = h pa ' ,pv . In fact, f\ can be 
obtained from f 2 by substituting h pvp ' a = h pvpa , and therefore is not independent. Another 
substitution of h pvp ' a = h pv > pCT + h PP]VCT + h vp ^ a into f 2 results in a new form, 

U = (PlPl) + (pIpD + (pIpD + 2(piP 2 ,PiP 2 ) + 2(piP2,pi +pl). (33) 
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Of course, this f± is not an independent form because this is derived from / 2 . But as will 
be shown below, this form is needed to express f^ a ■ For this reason, it is convenient to 
divide the whole space of the form f 2 into the form of f± and the rest, and regard f± as an 
independent form. By subtracting / 3 from / 4 , one obtains a simpler expression, 

h = (Pi.Pa) - (PiP2,PiP2)- (34) 

Thus there exist three independent forms, $2, fo, h, where f% denotes the rest of f% explained 
above. 



By substituting these forms into (IT7I) . one obtains 



K 



r equal / " \ J " J " f XJ " ^ 

J2R ~ U "77 ^ ~ OK-' 2 ' 



-i Hf =11 

57 9 /2 ~25" 
/ 3 e r=Q) (y/a - f/s) = §/a - |/ 5 , ( 35 ) 

requal I \ f f 

hR ~ \h) 3 /5 ~ 25 /5 ' 
Thus the eigenvalues of Rrk are 11/25 and 3/25. 
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